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On Certain Expansions of Elliptic, Hyperelliptic and 
Related Periodic Functions,* 

By F. R. Moulton. 



§ 1. Introduction. 

The functions treated in this paper are defined by differential equations 
of the types which arise in many physical problems. For this reason special 
interest attaches to the real solutions. The problem of finding the properties 
of the solutions and their expansion is, however, one of analytic functions, and 
therefore it is not convenient to limit the discussion to real values of the 
variables. Consequently the analytic character of the solutions is determined 
and, because of the possible practical applications, especial attention is given to 
developing convenient methods for constructing them. 

When the differential equations treated here are suitably specialized, they 
define certain elliptic or hyperelliptic functions. One of the most perfect theories 
in analysis is unquestionably that of elliptic functions. Nothing could be more 
systematic and elegant than the development of the properties of these functions 
when they are defined by the conditions that they shall be uniform and without 
essential singularities for all finite values of the independent variable, and that 
they shall have an algebraic addition theorem or that they shall be doubly 
periodic. These strong conditions make possible the detailed development of 
the theory, but they prevent its ready extension to more general cases. In this 
paper the point of view is altogether different. The functions are defined by 
differential equations, and the loss suffered in not explicitly prescribing their 
properties is partly balanced by the greater elasticity of the methods. 

In physical problems the questions of periodicity of motion and the asymp- 
totic approach of* the dependent variables to limiting values are of primary 
interest. In the case of elliptic functions there can not be more than two periods, 

* Bead before the American Matbematical Society under a different title, April 38, 1911. 
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only one of which can be real ; but in the case of hyperelliptic and more general 
abeliaii functions the number of real periods is not limited. In § 2 the question 
of the existence of the real and purely imaginary periods is treated. In this 
connection the qualitative character of the real curves defined by the differential 
equation is considered. In § 3 a number of special cases are treated, the simplest 
of which gives rise to the ordinary elliptic functions. In § 4 a method is given 
for expressing the periods by converging series in certain parameters of the 
differential equations. In § 5 the explicit construction of the periodic solutions 
from the differential equation is treated, and some of the important properties 
of the solutions are derived. In § 6 a method is given for constructing the 
solution from the first integral of the differential equations. In § 7 a method 
is given for constructing the solution, its properties being known from the 
discussion of §5, from the values of the dependent variable and its derivatives 
for special values of the independent variable. The whole problem of con- 
structing the solutions is there reduced simply to solving sets of ordinary 
linear algebraic equations. In § 8 application is made to the Legendre elliptic 
functions, and new expansions are given for them which, in § 9, are proved to 
converge for all finite real values of the independent variable if the absolute 
value of the modulus x^ is less than unity. The series are very convenient for 
practical use and are developed to the twelfth power of x inclusive. 

§ 2. The Differential Equations and Periods of the Solutions. 
Consider the differential equation 

where the J., are constants. Suppose the series converges if r<; \x\<C^B. This 
problem includes that of the motion of a particle moving subject to a^central force 
varying as any positive or negative power of the distance. The determination 
of the motion of a particle in the equatorial plane of an oblate or prolate spheroid, 
or in fact in the equatorial plane of any body whose surface can be expressed 
by zonal harmonics, depends upon an equation of the type of (l). Equation (1) 
also includes all those problems whose solutions are obtained from the inversion 
of elliptic and hyperelliptic integrals. Keference to any treatise on elliptic func- 
tions will show the abundance of these problems. 
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Equation (1) admits the integral 



x' 



= V/(x) = J ¥ Bjx^ + Blogx, (2) 

' ;; = -« 

where Bq is the constant of integration and where 

B=2A_i, Bj=^^^, y = — 00, , + 00 , except 0. 

The series /(as) also converges if r •< | a; | ■< 22. 

It follows from Cauchy's fundamental theorem on analytic differential 
equations that if, for any finite ti, the dependent variable x takes such a value ccj 
that the right member of (1) is regular, the solution of (1) is a regular function 
of t in the vicinity of t^ . Therefore the solution of (1) is regular for at least all 
finite values of t for which x stays in the ring r<C \x\<^B. It does not, of course, 
follow that for given initial conditions finite values of t exist such that x takes 
all values in the ring. 

Suppose the roots of / (a) = in the ring r <^ri< \x\ S^i-K ^ ^^^ 

tti, , ttp, where p is finite except in the trivial case /(a;) = 0. Suppose 

the a) are all distinct. Now join the circumferences r and i2 by a cut C which 
does not pass through an aj and which has no double point. From (2) we have 

/»* dx ,. 

It is well known from the general theory* that the values of t obtained by 
integrating from a; = cci back to a; = x^, oyer any path not passing through an o^ 
or crossing r^, Bi or C, are linear functions with integral coeflScients of Je — 1 
primitive periods Wg, . . . ., dp and /j, where Wg, . . . ., Op are the values of the 
integrals along closed paths around a^ and a^, aj and ag, . . ■ ■ , a^ and a^ 
respectively, and where I^ is the integral from to ai . These periods may not 
all be distinct and they are in general complex. 

We shall confine our attention now to the real and purely imaginary periods. 
Suppose the coefficients of (1) and (2) are all real. Suppose ai, . . . . , a^ are the 
real positive roots of f{x) = 0, and that /?i, • • • • , /?» *^6 ^^^ real negative roots 
of /(«) = satisfying the inequalities ri<\x\<Bi. It follows from the hypo- 
theses above that/(a3) is finite, continuous, and single-valued in the real intervals 
— Bi<:^x<^ — Ti, ri<;a;<^i2i, and it changes sign only at the aj and ^j. 

* Pulseux, Journal de Matli&matiquea, 1850, and Picard, TraiU cP Analyse, Vol. II, Chap. VIII. 
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Suppose / {xi) > 0, where a^ < JCj < a^+i , and that / (as) does not vanish for 
a« <C * <! <*i+i • Therefore x'l is real and we may suppose it is positive. The 
solution X will be a real function of t and x will increase with increasing real 
values of t until it reaches a^+j. It follows from (3) that this value will be 
attained for a finite value of t if a^+i is a simple root of f{x) = 0. If x = aj+j 
is a multiple root, x will approach a«+i asymptotically as t becomes infinite. 
Suppose a^+i is a simple root and that x has attained this value at ti^^. It 
can not pass beyond this value, for then f{x) would become negative and x' 
imaginary, whereas as is a regular function of t having real coeflBcients in the 
vicinity of this point. For the same reason it can not become complex as t takes 
increasing real values. It can not become a constant at this point, for d' is 

/ 




Fio. 1. 



distinct from zero, since a^+i is a simple root of / {x) = 0. Therefore as t 
increases beyond ^^^i, d becomes negative and x decreases continuously to a<, 
which it reaches in a finite time if a^ is a simple root of f{x) = 0. Then the 
sign of a;' changes and x increases to a^+i again. That is, under the conditions 
set down x oscillates between a^ and a^+i in the period 

r(<..,«.«)=2/;";^. (4) 

Similar results are true for all other closed intervals in the ring whose end-points 
are simple real roots of /(a!) = for which f{x) is positive. The periods are 
not necessarily distinct. 

In an interval in which /(a;) is negative we make the transformation 

t = a/ — 1 <r, 

after which the problem in t has the same properties it previously had in t 
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That is, there are purely imaginary periods corresponding to the closed intervals 
in which f{x) is negative, and for them x is real.* 

The solution involves two arbitrary constants, one of which is ^oi ^he other 
is additive to t, since (1) is unchanged by the transformation t=ztQ-\-t. The 
constant <o has no influence upon the analytic character of the solutions, but it is 
quite otherwise with Bf,. Since this constant is additive in f{x), changing it 
simply produces a vertical shift in the curves of Fig. 1. Suppose 5,, is increased. 
As a minimum approaches the jc-axis, the two real periods associated with the 
adjacent intervals for which fix) is positive approach infinity. On the other 
hand, the associated imaginary period approaches a finite determinate limit, as 
can be shown by forming the expression for the period and passing to the limit. 
If tti+i and tti+a are the two roots which approach equality, we have 

/(«)=—(«; — aj+i) (ai+2 — x)i> {x), 

where 4' («) > for a^+i < x < a^+j . Let the value at which a^+j and aj+g unite 
be a and -4' (a) = -4 > 0. Let a; = a + ^. When a^+a — a^+i is sufficiently small, 
•4- (a + ^) can be expanded in a power-series in ^ which converges for all x in 
a^+i < a; < aj+2 . Therefore the imaginary period is 

VAJa^,-a \/(a— ai+i)(ai+2— a) + (ai+i + ai+2 — 2a)^— P' 

where a = -7== . The limit of this integral as a^+i = a, at+z = a is ^=t . 

As ^0 increases beyond the value for which ai+i= a^+g, the purely imaginary 
period and one of the real periods giving real values of x disappear, and there 
remains in place of them a single period, which is real, giving real values of x. 
If the value of Bq decreases until a maximum of /(a;) crosses the a;-axis, the roles 
of the real and imaginary periods are interchanged. 

Suppose the positive roots of /(a!) = in the interval ri< |x|<jBi are 
tti, Og, as, oi^^ Suppose f{x) is positive in the intervals aittg, aga^ and negative 
in ag as . In the positive intervals x is real for real t, and in the negative inter- 
vals it is real for purely imaginary values o{t = V — It. In Fig. 2 these curves 
are given. Those which belong to real values of t are given in full lines and 

* While elliptic functions can have only two periods whose ratio can not be real, the functions defined by 
these more general equations can have any number of periods. This point has been insisted on by Casoratj 
in two memoirs in Ada MatJiemaiiea, Vol. VIII, 1886., 
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that which belongs to imaginary values of t is given in chain dots. In its case 
the ^axis should be read the t-axis. 

Since the solutions carry an arbitrary constant added to t, they may be 
displaced in the horizontal direction without changing their character. The 
figure has been drawn so that £c starts at aj at i = 0. It arrives at ag at i = J Ty^ 
and returns to aj at t=.T^. But if t takes purely imaginary increments when 
it arrives at |2^, x will describe the curve belonging to the interval a^a^. 
At if = i T^g -f- \/ — 1 J T^, X will have the value ag, and will return to ag at 




Fig. 2. 



< = i Tjg + V — 1 ^gj. But if t takes real increments after t = \T^'\- V — 1 J Tgg, 
X will describe the curve belonging to the interval a^ a^ , and will arrive at 04 at 

^ = 2 ^12 H" ^ — 1 2 ^8 H" 2 ^34 • 

The discussion of the branches for |£c| >-jB can not be made unless more 

properties of f{x) are given. 

§ 3. Some Special Cases. 

Case l. 5 = and only a finite numher of the B_j distinct from zero. 

If only a finite number of the B_j are distinct from zero, r = 0. If B is 
distinct from zero, the constant Bq is not defined for a? = 0. If a; has initially 
some value distinct from zero, and if for t = ti it becomes zero, ti is an essential 
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singular point. If B_j = for / = 1 and / >• m, and if B_jn ^ 0, the nature 
of the function in the vicinity, of < = <i, where x = for t = ti, can easily be 
determined from (2). We have at once 

TO 

x^ dx 



t — t 



1— »/n 



[I + Cix + ....]• 



TO +2 



(6) 



2 + 7nVi5_^ 
Inverting this series, we have 

« = P-^^t^5-TO(<-<:r]'^'{l + A(<-^i)'^+..-.}. (7) 

If the constants A^j are real, the coefficients Dj in (7) are real. The 
variable x considered as a function of t — ti has m + 2 determinations which 
vanish at t — i^ = 0. If «i is even, they are all complex, or two are real and m 
of them are complex for both positive and negative real values of t — t^ 
according as B_^ is negative or positive. If B_^ is negative, two branches 
are real for both positive and negative purely imaginary values of t — t^. 
If m is odd, there are m + 1 complex values of x and one real value for both 
positive and negative real values of t. In 'this case there is one real value of x 
for both positive and negative purely imaginary values of t. 

Suppose m is even, 5_„ is positive, and f{x) has simple zeros at 
a; = ai, a^, ^i, /^g, where ai^^a^-^H^B, — /^i-< — /8a<^. Suppose a5 = at 
f = ij . There are two real branches in the neighborhood of < = ^^ for both 
positive and negative real values of t — t^. Let us take t — ti positive and first 
consider the branch for which d becomes positive. It follows from (2) that on 
this branch a/ stays positive and x increases until it reaches the value ax . For 
real increasing values oit — t^, x' changes sign and x decreases to zero, which it 
reaches at a finite value of t — ti=- T{0, a^. At this value oit — t^, a;' changes 
sign by passing through infinity and x remains positive after passing through 
zero. Then the cycle to a^ and return is repeated. There is a similar cycle 
to /3i and return with the period ^(0, ^^ obtained by starting with the branch 
for which a/ became negative for real t — ti near zero. 

If t — ti takes purely imaginary increments after arriving at any odd 
multiple ot ^ T (0, %) on the branch which lies between and a^ , x will 
describe a real curve between ai and a^; and starting from the other branch 
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there is similarly a real curve between ^i and ^^ for purely imaginary incre- 
ments to t — <i after it arrives at any multiple, of T{Q, /3i). Fig. 3 shows the 
curves for real t — ti in full lines, and for imaginary t — t^ in chain dots. The 
curves are drawn satisfying the initial condition aj = at < = 0, and by con- 
tinuous variation of t from this value all the branches shown can be obtained. 

Now suppose that m is again even, that/(a;) has simple zeros at a3 = ai, ag, 
/3i, /?g, but that B_^ is negative. The only qualitative difference between 
this case and that just treated is that the full curves belonging to real t — t^, 




Fio. 8. 

and the chain-dot curves belonging to a purely imaginary t — ti, are inter- 
changed in position. 

Suppose m is odd, B^^ is again positive, and /(a?) has simple zeros at 
a; = ai , a^, ^i, ^2, where ai < ag < -B, — ^i<C — 1^2 <C + ^« ^^ t^^is case 
there is only one real value of x both when t — ti is positive and when it is 
negative, and in both cases x is positive. It t — t^ is a positive or negative 
pure imaginary, there is one real branch on which x is negative. For real 
t — ti, whether increasing or decreasing from t — ti = 0, x increases to ai and 
then returns to zero after a finite time T{0,a^. For purely imaginary t — t^, 
X decreases to ^i and returns to zero after the .finite time */ — 1 7(0, /Sj). If 
t — ti increases by purely imaginary increments from J ^^(0, a^), x remains 
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real and increases to ag, and then, decreasing, returns to aj at the finite time 
V' — 1 T{ai, otg). If t — ti increases by real increments from V — 1 J T(0, /Sj), 
X remains real and decreases to ^^ and then returns to ^i after the finite interval 
y(/8i, /Sg). Fig, 4 gives the real curves with the arbitrary determined so that 
35 = at < = 0. By continuous variation oft through real and purely imaginary 
values over the proper ranges, all the curves will be described. 

The case where B_^ is negative differs qualitatively only in that the curves 
for real and purely imaginary values of t — ti are interchanged in position. 




Case 2. B ■=■0 and only a finite number of the B^ distinct from zero. 

We now consider the case where jB = 0, B^^O and the -B^ = 0, y >• n. In 
this case iU = oo . The question of chief interest is the nature of a; as a function 
oft for those values of t for which a; = oo . Let jc = l/^ and (2) becomes 






r 



(8) 



If n < 4, the case of elliptic or exponential functions, ^ is a regular function of 
< at <i , where ^ = 0, a; = oo at < = <i . In other cases we find from 



that 



•^0 VB, + B,_,^+.... 
^ = [in^BAt-t,rf'{l^...,\ 



(9) 



(10) 



186 



MoTJLTON: On Certain Expansions of Elliptic, 



The discussion of the real branches of the curves defined by this equation is the 
same as that of (7). 

Suppose a is the largest positive root of /(cc) = and that /3 is the negative 
root having the largest modulus. Suppose further that these roots are simple. 
Then it follows from (9) that will vary from to 1/a or to 1//3 in a finite interval 
of time T{a, oo ) or T{^, cc ). Under the hypotheses which have been adopted, 
if t — ti increases continuously beyond T{a, ea) or r(/3, oo ), by real increments 




X 







Fig. 5. 

if Bn is positive and by purely imaginary increments if j?„ is negative, ^' will 
change sign and ^ will return to zero. For 5„ positive and n even the infinite 
branches of the real curves are given in Fig. 5. 

If 5„ were negative, the curves would be similar for purely imaginary 
values of t. 

If n is odd and 5„ positive, there are positive real branches for real values 
of t — ti, and negative real branches for purely imaginary values of t — t^. 
Qualitatively the figure for the real branches would be the same as Fig. 6, 
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except that the negative branches would belong to purely imaginary values 
of t — ti. If £n were negative, the results would diflfer qualitatively from those 
where jB„ is positive only in that the real and purely imaginary values of t — ti 
would be interchanged in position. 

Cask 3. / (as) a polynomial of the fourth or third degree. 

Suppose the roots of /(a) = are aj, ag, as, a^ and that the coefficient of a;* 

is unity. Suppose aj, , ai are real and distinct and that the notation has 

been so chosen that aj <; ag •< as •< at . We now make the linear fractional 

transformation 

^ « « « c 2 c c 

(11) 



85 — as a4 — a3_| — ^8 ^4 — ^a 



x — az a4 — as ^ — ^z'^i — ^a' 
which transforms ag, as, at into ^g, |s> ?t respectively, and any other point x 
forming a certain anharmonic ratio with Og, Os, at into a ^ forming the same 
anharmonic ratio with ^g, ^s, ^4- We assign the values ^2= — 1, ^8= + 1, 
^4=4- l/«, where « is determined by the condition that the point a^ shall be 
transformed into — 1/x. This is, of course, the Legendre normal form. Then 



(11) becomes 
Solving, we get 



(1 + xY _ ai — as at— ttg _ ^^^ 

To — . %y ^ 1. 

— xy «i — «2 *4 — "8 

The former of these two values of x is less than unity and will be adopted. 
It leaves the transformed points in the same order on the real axis as the 
original ones. Then equation (2) becomes 



^' = ^/^^/(l-P)(l-^2p)^ 
A — («a — «i)(«4 — «8) 



■ (13) 



{\-xY 

If /(») were of the third degree, that is, if at were infinite, the transformation 
would lead to the same result. 

In the case under consideration ^ is a regular function of t for all finite 
values of t. It has four periods, 

r(-l/«, -1), T(_l, +1), ^(+1, +1/;,), 7(1/^, +«,_oO,-l/«). 

The first and third are obviously equal. Since the transformation ^ = — l/xyi 
leaves the form of the first of (13) invariant, the second and fourth periods are 
24 
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also equal. One period is real and the other purely imaginary. The real 
curves are shown in Fig. 6. 

Case 4. The differential equation x' = ^/ {i — x^) {1 — x^ x% 

The case n >■ 2 will be discussed here only to show in what general respects 
the functions defined by this diflferential equation diflfer from those heretofore 

% 




Fig. 6. 

treated. For simplicity /(as) is taken in the Legendre normal form. But it 
is noted that if we should start with /(aj) having the zeros ai,a2,a8,a4, 
the transformation (11) would not in this case reduce the differential equation 
precisely to the normal form. 
The solutions of 

x' = >/{l—x^){\~x^x^) (14) 

are regular functions of t, certainly for all finite values, except possibly those 
for which »= d: 1, ;ea;= ±1, a; = ec 
points in turn. 



We shall examine the function at these 
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Suppose as = at < = 0. One of the determinations of the radical is real 
and positive. If we choose this one and take t real and positive, x will increase 
to unity, which it will reach at the finite time 

t, = T(0,l) = X'^==^—^. (15) 

If n were even and the negative determination of the radical had been taken, 
X would have arrived at — 1 at the time ti. 

We must investigate the nature of x as a function of t in the vicinity of 
t = ti. Let a; = l — ^; then (14) gives 

t — t -— r^ ^^ (\&\ 

''~ Jo ^^{2 — ^)il — ^+2xn — »^^')' 

Expanding the integrand as a series and integrating, we get 

-(<-<i) = ^'^[«o + ai^+ ••••], 
the Ut being real and a^ positive. Inverting this series, we have 

^ = 6iT + 6ar® + , 6i>0, 



(17) 



It follows from the second of these equations that r is not changed by multi- 
plying t — <i by any n-th root of unity. Therefore, there are n values of t — ti 
having the same modulus which give the same value of t. As t — ti passes 
through zero, f passes through zero; and its value after changing sign is, for a 

given modulus, equal to that before multiplied by e "~^ , which is a complex 
number except when «=2. Therefore, if f — ti is real any real curve belonging 
to it becomes complex when t — ^i passes through zero. It is obvious therefore 
that these functions can not represent the motions of physical bodies. 

If n is odd, | has two real values, one positive and one negative, for real 
negative values of t — t^. In this case ^ is not real for any other real or purely 
imaginary values oft — ti. If w is even, ^ has only one real value for real values 
of t — ti, and this value is positive and the same whether t — t^ is positive or 
negative with a given modulus. But in this case ^ is also real, though it is 
negative, for t — ti equal to a real positive or negative number multiplied by 
any n-th root of — 1. If w is 2, these are the purely imaginary values of t — ti 
encountered in the problems heretofore considered. In general, if n is even 
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and not congruent to zero mod 4, two of these values are conjugate pure im- 
aginaries. But if w is a multiple of 4, all these values of t — t^ are complex. 

It follows from this discussion that if we start with jc = 0, as' real and 
positive at < = 0, a; will take the value + 1 at a finite time < = <i . At this point 
there are n — 1 branches of the curves satisfying the differential equation, and 
there are values of t — <i, real or complex, for which x remains real and in- 
creases, and others for which it remains real and decreases. On the branch on 
which it increases it arrives at -fl/x at a finite time ^ = ^g. At this point the 
description of the branching differs from that at « = + 1 only in that positive 
and negative as applied to |, where a; = 1 — ^ as before, must be interchanged. 
Consequently there is a real increasing value of x beyond a; = 1/x for real 
positive values oi t — t^. If we follow this curve, x will tend toward infinity, 
which it will reach at the time 

^"/''-'=X. V(i-.'(i-«'.r ''*' 

which is infinite unless n <:^ 4. 

Suppose « = 3 and let x=. — 1/^. Then, if a; = oo at < = <8, we have 



'' = fo 



(19) 



the Ui being real and Uq positive. From this we get 
^ = boit-t,Y-^ 

»= ~ ^yi f ]co + Oi(< — <8)+ f, co>0. 

The discussions for — 1, — 1/x and « are essentially the same. We see 
that the details are quite different when n is even from those when it is odd, and 
that 85 = 00 for finite values of t only if n is 2 or 3. In the former they are 
regular points, and in the latter poles of order 3. For n = 3 the real curves 
are shown in Fig. 7. 

§4. Explicit Computation op the Periods. 
Suppose it is desired to find the period corresponding to the interval «», ai+i- 
One method is obviously to compute the numerical value of (4) by mechanical 
quadratures, as Legendre did in the case of elliptic integrals. But this is not 
sufficient in the present discussion, because in constructing the solutions in §§5-7 
it will be necessary to have the periods developed in converging series. 
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Since the zeros of /(») in ri<|{c|^-Ki are ai, . . . ., a^, we may write 

f(x) = (» — ai) . . . . (a; — a,,) -^ix), (20) 

where ^Csb) is finite and distinct from zero in the ring ri<|x|<^i. Suppose aj 
and tti+i are real and that /(a5o)^0 for ai^^o'C *i+i> ^o that the period is real. 
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Let 



Then (20) becomes 



Fig. 7. 



a 2 + 2 ^• 



/(»)=(i-anr^ 



+ «{+! — 2a^ , at+i — a< 



^];c(a 



(21) 



(22) 



where ;g(^) is finite and distinct from zero for all ^ corresponding to *'i<|«|^^i. 
Those o; which are complex occur in conjugate pairs. In these cases we take 
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their products and obtain factors of the forna ^j^^ + 7?? + ^^, where ^j, y, and 
hj are real constants. Then we may write (22) in the form 

/(«) = ^fe±=^'(l_a(l + px?) 

• • • • (1 + 9^m + U + MT . . . . [1 + 3\..^ + /«.P]i^(a, (23) 
where pi, • • • -, p*, Jli, i^i, • • • ., ^i, (ii and the coeflScients of F{^) are real con- 
stants, and where F{0) = 1. Moreover p^, • . • •, p^ are less than unity because 
tti and ttj+i are consecutive real roots of f(x) = 0. 

In order to compute (4) we must take the square root of the reciprocal 

of /(«). The square roots of the reciprocals of the factors (1 + pi|) (1 + pjl) 

are expansible as power-series in ^ which converge for all (^|Sl- The expan- 
sion of the square rooi of the reciprocal of 1 -|- Jl^^ + (ij^^ as a power-series in ^ 
converges if |^|-< 1/'V^|/%|. We shall suppose at present that the \(ij\ are all 
less than unity. Then (4) becomes 

where Pj, Pg, .... are polynomials in pi, , p^, X^, [ij and the coeflBcients 

of F{!^). The series under the integral sign converges for lll< 1. Since 

J_a ^T^Zp"*^ ~ "' J-i VT^:^^ ^ ~ 2^(«.' y ' 
the period is 

Since 

(2n).^ _ 1 .3.5 2n— 1 ^ 

22«(„/)2— 2.4,6 2n ^ 

and the integrand of (24) converges for |^| •< 1, the right member of (26) is a 
convergent series. 

In the simple elliptic case 

/(a!)=(l-jB^)(l-x«A x^<:i, 
the real period is 

r=^{i + (l^)V + .... + (i^i^^^^g=i)V + ....}. (26) 

Suppose now that \(ij\ > 1, but that all other I/KjK 1, arid that 
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where /3 and y are real. The integral in this case can be computed from — 1 
to ^, and then from /3 to + 1. Let 

Then >7 = — 1 when | = — 1, and >; = + 1 when | = /3. The quadratic 

(1 4- 8Y 
factor [(^ — ^y + y2] becomes ^-— J-^ (1 — )7)^ + y^- The square root of its 

reciprocal is expansible as a converging power-series in >? if 1 17 1 <[-J 1 -|- ^ _ ^ . 
Therefore the first of the two integrals is 

y/„^^ ^^^ r+' ll + Qiyi+ Qzyi'+'-'-ldy; 

J-1 ^l—yj 

For the second integral there is a similar transformation and expansion. 

If two of the I iij I were greater than unity, the integral in x from a^ to a^+i 
would be broken up in a similar way into the sum of three integrals which 
would be defined by converging series. This process of division can obviously 
be extended to include the case where any finite number of |/m^| are equal to 
or greater than unity. 

§5. Explicit Constritotion of the Periodic Solution. 

Suppose it is desired to find the solution for which a; oscillates between 
ttj and tti+j. By the transformation (21) the diflferential equation becomes 

^'=4V(i-p)(i+p,^). . . .(i+p.^)[i+;ii^+f*iP] . • • • [i+^iRM^^ll). (27) 

where ^(^) is a power-series in ^ which converges if |^|<1. 
For convenience let us make the transformation 






(28) 



where the Cj are coefiBcients of expansion of F{^). The parameter (i will guide 
the computations ; in numerical applications it will be taken equal to unity and 
the bars on the p^, Jlj, /Uj and c^ will be omitted. After the transformations (28), 
equations (27) and (25) become 



f = 4^/(l-P)P(^J|t^), 



(29) 
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where P is a power-series in (i which converges if |^|<1, provided all |/M;|-< 1, 
and where Q{(i) is the series (25). The coeflScient of fi^ in P contains ^* as a 
factor. 

In order to place in evidence the period in the construction of the solution 
let us make the transformation 



< — <o = 3 ©(/«)*• 



(30) 



Therefore, the period in the variable r is 'in. If we denote by dots the 
derivatives with respect to r, we have 



dP 



(31) 



^ = e(^)V(l-P)P(^;/«), 

By Poincar^'s extension of Cauchy's theorem* the second of these equations 
can be integrated as a power-series in [i, and | ft | can be taken so small that 
the series will converge for any preassigned interval for t. The same result 
is eflFected by taking p^, \, ^ jand q small, and (i equal to unity. 

Let us construct the actual solution of the second of (31) with the initial 
conditions !^\0) = — 1, ^ (0) = 0. We have 

P =l + Pi^^ + P,P(«« + 

dp 



d^ 



Pi/t -f 2P2?/t*=' +.... + « p„r-v" + •• 



(32) 



where the P„ and Q^ are known constants, and the ^„ are functions of r to be 
determined. Substituting (32) in (31) and equating to zero the coeflficients of 
various powers of fi, we find 

lo + lo = 0, 



Zn "^ Zn — ^n(?0> • • • • > Sn-l)> 
> 

where Rn is a known polynomial in |o, , In-i • 



(33) 



* Let Methodes Nouvelles de la M6eanique GiUste, Vol. I, p. 58. 
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Equations (33) can be integrated in the order in which they are written. 
Since the initial conditions are identities in fi, we have 

With these initial conditions the solution of the first equation of (33) is 

^0 = — cos t. (35) 

Then the second equation of (33) becomes 

li + ^i = -iPi-iPiooB2t, 

whose solution satisfying (34) is 

^i = -iPi+iPi0O82'e. (36) 

Let us suppose that ^o , , ^^-i have been computed and that it has been 

found that |^, y = 0, , n — 1, is a sum of cosines of multiples of t, the 

highest multiple being y+1, and then consider the right member of the (n + l)-st 
equation of (33). Its general term, aside from a constant coeflBcient, has the form 

[K]=eK.^•••^^ (37) 

It follows from (31) and (32) that the exponents and subscripts satisfy the relations 

h + iz+ + *n = 45« + 1| 

hJi + hh+ '••• + K3n <n — (4_i). 
Bach element of the general term being a sum of cosines, the product is a sum 
of cosines. The highest multiple is 

Nn = ilUl + 1) + hUz + 1) + • • • • + inUn + O, 

which, by (38), has as its greatest value n + !• Therefore, the differential 
equation at the (n + l)-st step is 

'L + L = A""^ + ^f ^ cos -r + .4(") COS 2 -r + . . . . + ^JTii cos (» + 1) t. (39) 
The solution of (31) is periodic with the period 2 7i in t for all \(i\ 
sufficiently small. Therefore, 

where the identity holds in both fi and t. It follows from the fact that it is an 

identity in (i that 

^j{t+2n) = ^,{r),j = 0, ...., «. (40) 

Since this relation is an identity in t, each ^, is separately periodic with the 
period 2n. 
25 



(38) 
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Now consider the integration of (39). Since the solution is periodic, the 
coeflScient J.^"^ is zero. This condition has been secured at each step by the factor 
Q^{{i) in the second of (31). With A!f> — 0, the solution of (39) satisfying (34) is 

^„ = a^"^ + af^ cos T + af > cos 2 r + + ajfjj cos (n + 1) t, 



ar = ^i"', 



A(n) 



af> = — ^^»> + 2 



n+l Mn) 



(41) 



y=2 j"- 1 

These results have the properties assumed in the beginning, and the induction 
is therefore complete. 

§ 6. Direct Construction of thk Solution prom the Integral. 

The properties of the solution being known, it can conveniently be found 
from the first equation of (31). Let this equation be written in the form 

e-Q'{^){i-e)P{^)^)=F{i^)(i)=o. (42) 

This equation holds identically in ^i and t. Since and | are expansible as con- 
verging power-series in j«, F can be expanded as a power-series in ^ of the form 

i?' = i?'o + i^i/it +....+ i^„(«" +....= 0. 

Since this equation is an identity in /«, we have 

-P'n(ii,6) = 0, y=0, ...., n; « = 0, ....,00. (43) 

Now suppose the expressions for ^j and |^ given in the first of (41) are 
substituted in (43). It follows from (42) that the largest value of j in F^ 
is j = n. The result of the substitution will be a sum of cosines of integral 
multiples of r, the highest being n + 2. Therefore, (43) may be written 

F^ = 5^">+ -Bf) cos -r + .... + B'^^ cos (n + 2>r = 0. (44) 

Since this equation is an identity in r, we have 

5j«>=0, y = 0, ....,nH-2. (45) 

It follows from the form of (42) and the fact that in P(^, (i) every term except 
the first, which is unity, is multiplied by (i that the 5j"' have the form 

5^") = — 2ai"> + Gf>, 

Bf = — 24«> — 3c4»' + Gf\ 

^in) =^ _!_ 4ai"^ + <7'*' 

5f) = + (/- 2)ai!i\ - (i + 2)ajf, + C7f> , /= 2, . . . ., n, 
mi= + (n--^y^' + + 0^1, 



(46) 



Hyper elliptic and Related Periodic Functions. 197 

where the (7j"*, ^ = 0, , n + 2, are known functions of af\ . . . ., af~". 

Therefore if af\ , aj"~" have been computed, equations (45) determine 

(4"', ••..,aJi'^i. The solution is most conveniently made by starting with the 
last equation and using them in order to the first, with the exception of the third, 
which is redundant. If the constant of integration in (42) had not been deter- 
mined on the start so as to belong to the solution with the period 2n, the 
coeflBcients of the (n + l)-8t term of its expansion as a power-series in (i would 
be determined by the third equation of (46). 

§ 7. Direct Construction op the Solution from the Values of ^ and its 

Derivatives at t = and r = n. 

We may write (31) and their successive derivatives in the form 

^ = i2,(^, (i), 



where the Ej{^, (x) are power-series in ^ and fi. The derivatives obtained in 
the right member are of course eliminated at the successive steps by the first 
equation. Since f= — 1 at'r = and ^=-f-l atr=:7i identically in fi, 
we have 

^(o) = i^,(oV=-i, 

^=0 



}=0 



(48) 



^(0) = 2 ^Xoy = i?n(- 1, i«) = 2p<f(- l)(i^, n = 1, . . . ., 00 

l{n) = i 'Un)f*> = Rn{+ 1, /«) = ^m+ ly. 

Now consider the determination of the coeflScients of ^j. Its form is given 
in (41), from which it follows that all odd derivatives vanish at both r = and 
r = »t. All Rzk+ii^t l^) *l80 vanish at t:=0 and r=zn, because they carry 
v^l — ^2 as a factor. Hence the coeflBcients of equations (48) are identically 



(49) 
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zero when n is odd. Since equations (48) hold identically in /*, we have for 
even values of n 

^^(o) = o, y=i, ....,00, 

|,(7t)=0, 

2n 

6(o)=i>ii'(-i), y=o, ....,00, 

2n 

or, making use of (31), 

a\,^ -\-af + .... + aj$i = 0, 

(4^ - ai^> + . . . . + (- l)^+'aj4'i = 0, 

+ ai^ + 2^"4-'' + . . . . +(y+l)^"aji>i = i)^(- 1), «=1, 2, . . . . , 00 , 

- aP + 2^4^ _....+(_ i)W (y + i)2n„j^^ =i>^^(+ 1). 

Since the right members of these equations are known, the a'i' are uniquely 
determined by this infinite set of equations, provided the determinant of the 
coeflBcients of any y + 2 of them is distinct from zero. It is clearly necessary 
to use at least one of the first two equations, which alone involve aj*. If those 
equations are taken for which t = 0, the determinant is an alternant and is dis- 
tinct from zero. Therefore the solutions are uniquely defined by the properties 
which have been established and by these equations. This is the simplest 
method of computing them. The results of the computation are easily checked 
by substituting them in any of the infinitely many superfluous equations. 

§8. Application to the Legendbe Elliptic Functions. 

The results which have been obtained in §§ 5-7 include the Legendre Elliptic 
Functions as simple special cases. Consider sn r, the sine amplitude function, 
which is defined by 

x" = _ (1 + ^)x + 2xV. J ^ ' 

The parameter x^ plays the role oi ft above, and the solutions are in powers of x^. 
In order that the results may agree with the usual notations, we take the initial 
condition a5(0) = 0. In this case we have 

g(^)=g(x^)=n-(i)v+ (i^^)V + .... + Q;^;;;;^;;~7 ^^"+ (51) 

Then 

»= g(x2)x/(l-a^)(l— x^a^). 
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The solution of this equation as a power-series in x^ is 

a; = Kq + JBgX^ + a;4X*+ (52) 

Making the computation by the methods of §§ 5-7, it is found that 

Xq = sin T, 

X2 = tV [sin T + sin St], 

Xi = TS4 [7 sin T + 8 sin Sir + sin 5t], 

ajg = 1V3 [67 sin T + 82 sin St + 16 sin br + sin 7^], 

ajg = tV* [738 sin t + 945 sin 3t + 230 sin 5t + 24 sin 7t + sin 9r], l (53) 

Xio=^Ti, [8808 sin t + 11,661 sin S-r + 3264 sin 5t + 442 sin 7r 

+ 32 sin 97 + sin ll<r], 

a;i3= yV« [110,728 sin T + 150,486 sin S-r + 46,519 sin 57 

+ 7440 sin 7t + 718 sin 97 + 40 sin 117 + sin 137], 

Bach of these expressions carries sin 7 as a factor. Since x(n/2) = 1 for 

all 0^, it follows that x^{n/2) = 0, n=l, , 00. Therefore each x^Jn/i), 

« = 1, . . . . , 00 , carries cos 7 as a factor. Since sc(7t/2) also vanishes for x^ <[ 1, 
it follows that each x^^ ^^^ co8^7 as a factor. After sin 7 cos^7 is removed, the 
result can conveniently be expressed as a sum of even powers of cosines of 7. 
The explicit results are 

Xq = sin 7, 

ajg = I sin 7 cos* 7, 

Xi = TJlr sin 7 cos* 7 [5 + 4 cos* 7], 

jcg = -g^ sin 7 cos* 7 [10 + 11 cos* 7 + 4 cos* 7], 

ajg = t1^ sin 7 cos* 7 [389 + 500 cos* 7+ 272 cos* 7+ 64 cos^ 7], 1 (54) 

a;j„= if^T'sin 7 cos* 7 [1054 + 1499 cos* 7 + 984 cos* 7 + 368 cos® 7 

+ 64 cos* 7], 
jcjgzz TTrVT'iSin7 cos* 7 [12,242 + 18,759 cos* 7+ 13,800 cos*7 

+ 6448 cos«7 + 1856 cos«7 + 256 cos" 7], 



These expressions are very convenient for use in practical computations. 
Let y = en 7, the cosine amplitude function. Then we have 

[2»o + ia!8^«*^]ia-8yX*^' 



y = V 1 — a* = V 1 — ag J 1 



1 — xl 



(55) 
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It follows from this equation that y is expansible as a power-series in x^ and 
that the coefficient of x*" is a sum of cosines of odd multiples of r, the highest 
being 2n + 1. The explicit expression for y can be derived from this equation, 
but it is just as conveniently found from the dififerential equation which y 
satisfies, viz., 

y-q{x^) V(l-y2)(l-x2 + x2y2), (56) 

where y(0)=l, y(n/2) = 0. Forming the successive derivatives, equations 
corresponding to (49) are obtained. Since y(7t/2) = 0, it follows that each 

yg„, n = 0, , 00, carries cost as a factor, and since y(0)=l, y (0) = 0, 

it follows that each y8„, n = 1, . . . . , oo , carries sin^r as a factor. These results 
also follow directly from the properties of x and the form of (55). Making 
the computation, we find 

y = yo + yzx^ + ytx* + — , 

yo = cos T, 

^2 = — T COST sin^T, 



^4 = — AcosTsin^T[7 + 4cos^t], 

2/, = — ifiTCOST sin^T [17 + 13 cos^T + 4 cos*t], 

2/8 = — tIts COST sin* T [769 -f- 660co8*T+ 304co8*T+ 64co8*t], 

2^10= — ifiT^cosTsin2T[2289 + 2149cos*T-f 1192co8*T + 400co8*T 

-4- 64 cos^t], 

yi2= — TirVj» COST sin* T [29,023 + 28,727 cos^T + 17,800 cos*T 

+ 7472 C08«T + 1984 cos«T + 256 cos^"t], 



(57) 



Now let 2 = ^1 — 
satisfies the equation 



x*a;'= 



the delta amplitude function. The variable z 

(58) 



a = Q (x^) V (1 — 2*) (— 1 + X* + 2*). 
Making the computations as before, we find for the explicit value of z 

Z = 1 + ZgX^-f 24X* '+ , 

Z2 = — Jsin^T, 

Z4 = — I sin^T [1 + cos*t], 

Zg = — TiVsin^T[4 + 5cos^T + 2cos*t], 

Zg = — jVs sin*T [40 + 56 cos^T + 32 cos*T + 8 cos't], 

Zio=: — TTso sin^T [448 + 677 cos*T + 464 cos*T + 176 cos't 

+ 32cos«t], 
Zi2= — ji^ii sin* T [1344 + 2152cos*T+ 1586 cos*T-f 756 cos't 

-h 224 cos*T + 32 cos^°t]. 



[ (59) 
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The functions sn r, en r and dn r arise very naturally in many problems of 
mechanics and geometry.* The series (54), (57) and (59) are convenient for use 
in solving these problems, and the coefficients which are given here can be relied 
upon as being correct. But the point of chief interest is that if the diflferential 
equations defining the functions are of the more general type considered here, 
the same methods of obtaining their explicit values can be used. 

§9. The Realm of Convergence of the Series (52), (57) and (59). 

For simplicity let us first consider the series for the sine amplitude function, 
which is defined by 

x = Q {x') '>^ {1 — x^) (1 — x^ x^), X = at r = 0, 
x = — Q^ (x2) {l+x^)x + 2 Q^ (x'') x^x 

It follows from the second of these equations and the fact that Q (x^) converges 
for all I x^ I <| 1 that the solution for x considered as a function of x^ is regular 
for all I x^ I <H and all finite values of t, provided x does not become infinite 
for some value of .t. In the solutions (52) and (54) r enters as a parameter and 
takes only real values. Therefore, if x does not become infinite for any real r 
on the interval to 27t, the solution (52), being periodic with the period 2n, 
converges for all | x* | <[ 1 and all finite real values of t. 

In order to establish the result that x does not become infinite for a real 
value of T, we shall show that the integral 

dx 



(«») 






is not real. It is necessary to specify the path of integration from to « . We 
shall integrate from zero to 1 — e along the real axis; then along the circle with 
radius e and center at 1 to the line joining 1 with the point 1/x, which will in 
general be complex ; then along the line from 1 to 1/x to a distance e from Ijx ; 
then along the small circle e around Ijx to the line from Ijx to oo whose incli- 
nation to the axis of reals is the argument of 1/x; then along this line to oo . 
This may be indicated by 

j=j,+/,+/8+/t+/5=/+ J + /'+ r +r, (61) 

where /j, . . . ., /s represent the integrals (61) in order. It is easy to show that 
^^ L = ^"^ Z = 0. The transformation x = Ijx^ throws L into the form L. 

* See Appell and Laconr, Fonctioni Elliptiquet, Obap. V. 
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Therefore we have 

j_ „ r^ dx r^^' dx . . 

~ -^0 g(x2) V (1 — »^) (1 — x'x') -^1 Qix") >/{l — a?) {1 — x'ai'y ^ ^ 

By the definition of Q{j^) we have 

f* dx n Q{x') _ n 

Jo g(x2)V(l— a;'')(l — xV) 2 Q(x^)~ 2 ' 



which is real whatever the value of x^. It is one-fourth of one period of the 
elliptic function. The other integral is one-half of the other period, and there- 
fore it can not be real, since the ratio of the two periods can not be real. 

Now consider any other path of integration to infinity. If none of the 
branch-points ±1, ± 1/x is included between the two paths of integration, 
the value of the integral is the same. Suppose only the branch-point +1 is 
included between the two paths. Then the new path can be deformed, without 
crossing a branch-point, into the old one plus an integral from zero to + 1 and 
back, after making the circuit of this point. Hence, in this case the value of the 
integral is 4J1-I-/3, which is not real. Suppose only the branch-point 1/x is 
included between the two paths. The new path can be deformed, without 
crossing a branch-point, into the old one plus an integral from 00 to l/x and 
back, after a circuit around 1/x, along a line whose inclination to the axis of 
reals is the argument of 1/x. Hence, in this case the value of the integral is 
2/1 4- -4 — 2/1, which is not real. It is easy to see that the integral from zero 
to infinity along any path is 21^ + 1^+ 2mili + irrisla, where »Ji and wig are 
positive or negative integers, or zero. This sum can not be real and it follows 
that X can not become infinite for any real finite t. There/ore, the serieb (52) 
converges for all real finite values of 1, provided | x^| <; 1. 

An analogous discussion proves the same result for the expansions of the 
cosine amplitude and the delta amplitude functions. 

The Univbbsitt or Chicago, April 27, 1911. 



